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Abstract 

In this paper we have considered the problem of parametric sound generation in an acoustic 
resonator filled with a fluid, taking explicitely into account the influence of the nonlinearly generated 
second harmonic. A simple model is presented, and its stationary solutions obtained. The main 
feature of these solutions is the appearance of bistable states of the fundamental field resulting from 
the coupling to the second harmonic. An experimental setup was designed to check the predictions 
of the theory. The results are consistent with the predicted values for the mode amplitudes and 
parametric thresholds. At higher driving values a self-modulation of the amplitudes is observed. 
We identify this phenomenon with a secondary instability previously reported in the frame of the 
theoretical model. 
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I. INTRODUCTION 



Ultrasonic resonators are devices that confine the acoustic fields in a finite region of the 
space. When driven by an external energy source (e.g. vibrations of one of their boundaries), 
the reflection and resonance conditions imposed by the boundaries provide the conditions to 
obtain high amplitude fields, favoring the development of nonlinear effects, as the emergence 
of frequencies different from that of the driving. 

The knowledge of the field evolution in resonators in nonlinear regime is important both 
in its fundamental and applied aspects. One of the simplest configurations consists in two 
plane and parallel walls (an acoustic interferometer) where nearly one- dimensional standing 
waves are formed along the cavity axis. The system was the basis, forty years ago, of 
the first experimental observation of acoustic parametric excitation in a fluid The 
phenomenon, first observed by Faraday J4] and later described by Lord Rayleigh, consists 
in the emergence of oscillation modes at frequencies smaller than that of the driving, when 
a parameter of the system is varied periodically in time. The phenomenon is universal, and 
has been demonstrated in a variety of physical systems In the case of the acoustical 
interferometer, the length of the cavity (and thus the eigenfrequencies of the normal modes) 
is the time-dependent parameter, and the parametric excitation is achieved when the input 
energy is high enough to overcome the dissipative losses. The parametric fields usually 
appear as doublets, whose frequencies fi and / 2 add to match the driving frequency / , i.e. 
fo = fi + /2, although in some circumstances the half- frequency mode is observed. 

A theoretical description, based on the Mathieu equation, has been successfully applied 
to the description of these processes in an acoustical interferometer, allowing to predict the 
subharmonic spectrum and its excitation threshold [5, 1(5]. This approach has however a 
restricted validity, since it can not describe the further evolution of the parametric fields 
above the threshold and, on the other hand, it ignores the unavoidable effects of the higher 
harmonics of the driving, which are the main signature of non- or weakly-dispersive acoustic 
systems. 

The spectrum of higher harmonics can be controlled using additional dispersion mecha- 
nisms, such as bubbly fluids or walls with selective (frequency dependent) absorption 



Ostrovsky et al 



lOj ] designed a waveguide cavity where dispersion is introduced by the lat- 



eral boundaries, and the higher harmonics could be completely inhibited. As a theoretical 
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framework, a set of evolution equations for the trial of interacting modes was considered, 



resulting in an excellent agreement with the measured amplitudes a 
This theory has been also applied to the interferometer case 



Dove threshold. 
HQ , but the agree- 



ment with the experiment [7| was mainly qualitative. The discrepancies can be interpreted 
in terms of the influence of the first higher harmonics (those with higher amplitudes) on the 
parametric process, introducing additional features that can not be captured by the fully 
dispersive model. 

In this work we present a theoretical description of parametric sound generation in a 
fluid-filled interferometer, extending the previous models to include the coupling with the 
second harmonic. The stationary solutions both above and below the threshold are obtained 
and discussed. The solutions are compared with experimentally measured values, with an 
excellent agreement. The novel effects induced by the second harmonic (e.g. bistability 
or hysteresis) are discussed. Finally, we also report the existence of low frequency self- 
oscillations, occurring at pump values beyond the parametric instability threshold. This 
phenomenon is discussed in terms of secondary instabilities of the proposed model. 

II. THEORY 

The acoustical interferometer considered in this paper is composed by two parallel and 
solid walls, with thicknesses D, located at a distance L from each other, containing a fluid 
medium inside. Each medium involved in the model is acoustically characterized by its den- 
sity p, bulk modulus n, and sound velocity c, related as c 2 = -. The resonance modes (eigen- 
frequencies) of the resonator depend on these parameters, defining the acoustical impedances 
z = pc of each section. In the ideal (loseless) case, corresponding to an infinite impedance 
of the walls, the resonance modes are equidistant and obey the equation 

tan(&/L) = (1) 

where kf = 2nf /cf. However, in a real system the impedances have finite values, and 
the spectrum of the resonator is no longer equidistant, but distributed according to the 
transcendental equation 

27Z tan {k w H) 

where k w = 2nf /c^and 71 is the quotient between wall and fluid impedances. 



3 



Now consider the driven system, assuming that one of the walls vibrate with frequency 
/o- Then, above the parametric generation threshold the spectrum inside the resonator can 
be decomposed in two sets: the subharmonics resulting from the parametric instability and 
the higher harmonics of the driving, n/o, with n an integer number, resulting from the 
weak dispersion in the system. The amplitude of any of these harmonics decrease with the 
detuning, defined as the difference between the corresponding field frequency u n and the 
frequency of the closest cavity mode, oj c n , i.e. 5 n = u) n — u) c n . The energy transfer into a 
harmonic is then more effective under resonance conditions. In an ideal cavity, a resonant 
driving implies that all higher harmonics are also resonant with a cavity mode, enhancing 
the cascade of energy from the driving into many higher frequency components. However, 
an initial detuning 5q implies that the detuning of the higher order modes increases linearly 
with the frequency (Fig. 1(a)). As a consequence, the amplitude of these modes is reduced 
with respect to the resonant case. This situation can be more dramatic in a real system, 
where the mode distribution is non-equidistant and the detuning of the second and higher 
harmonics can be tuned to be much larger than that of the driving and the subharmonics, 
as shown in (Fig. 1 (b)). 

The above arguments suggest that a theoretical approach to parametric sound generation 
must take into account the effect of the second harmonic, whose amplitude is not negligible 
in general, but one can neglect the influence of the third and higher harmonics, assuming 
that they are sufficiently out of resonance. To simplify the analysis, we consider that the 
parametric generation is degenerated, i.e., the half harmonic of the driving is excited. 

Under these assumptions, and following the technique described in detail in [13] for the 
fully dispersive case (neglecting higher harmonics), the following system of equations for the 
evolution of the slowly-varying amplitudes for each mode is obtained: 

dP , . , x „ . „ / „o „ „^ c 



( 7o + i6 )P - if3 (P? + P 2 P*) + -P 
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dt wu u/ u ' v 1 ' L 
^- = -{lx + ih)Pi-i^P*Pl (3) 

d P 

^ = -(7 2 + ^ 2 )P 2 -^P 2 , 

corresponding to the fundamental (driving), subharmonic and second harmonic respectively. 
Other parameters are: P in the driving amplitude (pump), r y i the decay rates of each mode 
in the cavity, f3 = (1 + B/A) is related to the nonlinearity parameter B/A of the fluid 
and 5i = uji — w\ are the detunings. 



The dynamical system given by Eqs. ([3]) can be reduced to a simpler, dimensionless form, 
defining the new variables A = i^Po, A x = i ^ 2 ^ oll pi, A 2 = ij^Pz, E = i 2L ^° 71 p in , and the 
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parameters 7 = and A = S^/j^ With these changes we obtain 



7 7 2 

To ^ = " (1 + *A ) A + E - V - A 2 A*, 

1 - 1 1 ^ = -{l+i& 1 )A 1 + A Q A{ ) (4) 
7 2 ^ = - (1 + *A 2 ) A + ^A 2 . 

Equations (03) admit two different stationary solutions. When the pump amplitude E is 
below the parametric threshold, one can set Ax — 0. Neglecting the temporal derivatives in 
Eqs. 01]) we obtain 



2 r' 



^ | (l+A 3 + ^°i (7|Aol ++ J(l-A A 2) ) (5) 



f + A> 



\A 2 \=^-^\A 



1 1 2 1 1 2 

where we have defined Jo — I A) I an d -^2 = |A I • 

The solution given in Eq.flH]) reflects the fact that the amplitudes of both the fundamental 
and the second harmonic grow with the pump amplitude, as expected. The most interesting 
feature is however the emergence of multivalued solutions, i. e. the system shows a bistable 
or hysteretic behavior, when some conditions on the parameters are fulfilled. The condition 
for multivaluedness is found by imposing the existence of an inflection point in the S-shaped 
curve given by Eqs. ((Sj). This occurs when = 0, i.e. at 

IV = JJlz^M. ( 6 ) 

Note that the existence of bistable solutions requires that A0A2 > 1, i.e. both detunings 
must have the same sign and exceed a certain critical value. Figure 2 shows the monos- 
table (a), critical (b) and bistable (c) cases of the fundamental mode, for different sets of 
parameters. 

At a pump value given by 



E th = J(l + A?) [(1 + 7) 2 + (A - 7A 2 ) 2 ] (7) 
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the solution given by Eq.Q becomes unstable, denoting the threshold of the parametric 
instability. The new solution corresponds to a non-zero value of the subharmonic amplitude, 
given by 



l^l 2 = (-1 + AoAO -7(1 + AiA 2 ) ± y'p 2 - [A + A 1 + 7 (A, - A 2 )] 2 (8) 
1+A 2 

where we have defined 7 = 7 1+A ^ . 

Above the instability threshold, the fundamental mode saturates to a constant value 

|A | 2 = 1 + A 2 , (9) 

and A 2 is given by Eq.([5]). 

Note that in the limiting cases of large second harmonic detuning A 2 or large losses 7 2 (in 
both cases 7 tends to zero), the solutions (J5j) -(HI reduce to those obtained in previous works 
(see e.g. [131 ]) for the fully dispersive cavity. Physically a mode with large detuning or losses 
is strongly damped in the cavity, so in practice it can be neglected from the very beginning. 

The parametric instability can be either supercritical or subcritical, resulting in monos- 



table or bistable subharmonic respectively. This fact, different to the bistabi 
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ity of the 
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fundamental mode, is also present in the fully dispersive cavity, as shown in 

Finally we note that, in order to observe the bistable regime of the fundamental mode, 
it must occur with precedence to the parametric instability. Combining Eqs.([6]) and (J9j) , we 
find that the fundamental mode presents bistability whenever 

< |- (A A 2 - 1) - A 2 < 1. 
37 

III. EXPERIMENT 



A. Description 

The resonator consists in two piezoceramic discs (p = 7.70 x 10 3 kg/m 3 , c =m/s) with 
radius 1,5 cm, and thickness 1mm and 2 mm (corresponding to resonance frequencies around 
2 MHz and 1MHz, respectively), mounted in a Plexiglas tank containing distilled and de- 
gassed water (p = 1.00 x 10 3 kg/m 3 , c = 1480m/s at T = 20°C). Both sides are located 
at a variable distance L, and its parallelism can be carefully adjusted to get a high-Q in- 
terferometer. One of the piezoceramics, that with resonance frequency around 2 MHz is 
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driven by the signal provided by a function generator (Agilent 33220) and a broadband RF 
power amplifier ENI 240L. The experimental setup is completed by a needle hydrophone 
(TNUA200 NTR Systems) to measure the intracavity pressure field. The Fourier transform 
of the acquired signal allows to determine the spectral content of the resonator, and from it 
to quantify the amplitude of any frequency component. 

In this way, by changing the amplitude of the driving source, we are able to follow the 
bifurcation diagram of the resonator for a given set of parameters. Although the pump 
value and the decay rates r ) i can be unambiguously determined (the latter by measuring the 
line width of the cavity modes), the detuning parameters Si are difficult to control, and in 
general vary from one set of measurements to the next. The reason is the dependence of 
the cavity resonances with the temperature of the medium, which changes with time due to 
the external sources (heating of the driving transducer or ambient temperature variations). 
However, we have monitored the temperature variations to ensure that the detunings did 
not changed appreciably withing one set of measurements, so we could apply the theory of 
the previous section to the description of this problem. 

B. Bifurcation diagrams 

In the series of experiments we have obtained different bifurcation diagrams. They differ 
in the spectral content (frequency of the subharmonic pairs, weight of the higher harmon- 
ics) corresponding to different sets of detunings, but also in the behavior of the individual 
amplitudes. We illustrate two of such cases in Figs. 3 and 4. 

Figure 3 was obtained for a driving frequency fo = 1.991 MHz and L = 3.5 cm. In this 
case a nearly-degenerated subharmonic pair, with frequencies f\ = 1.018 MHz and fi = 0.974 
MHz, is excited. As shown in Fig. 3(a), corresponding to a driving amplitude of 25 V, the 
pair is highly asymmetric, the amplitude of one of the components, f2, being negligible with 
respect to f ± . The experimental bifurcation diagram is shown in Fig. 3(b) with symbols. 
The first evident feature of the diagram is the non-linear growth of the fundamental mode 
amplitude as the driving amplitude V in is increased, denoting a non-trivial influence of the 
second harmonic. At Vi n = 24 V, the threshold of parametric generation is reached, and 
the amplitude of the subharmonics begins to grow. Note that in Fig. 3(b) this amplitude is 
magnified five times for a better clarity. 
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The solid lines represent the analytical stationary solutions given by Eqs.flH]), © and 
The theoretical amplitude and pump values have been scaled by respective numerical factors 
taking into account the efficiency of the transducer, the sensitivity of the hydrophone, and the 
additional normalizations leading to Eqs.fllJ). We used the medium parameters p = 1.00 x 10 3 
kg/m 3 , c = 1480 m/s and a = 0.875. The decay rates were obtained by measuring the line 
width of the modes, resulting 7^ ~ 5000 s _1 . Finally, the detunings have been chosen to 
get the best fit to the experimental data. In Fig. 3(b) we used A = —2.35, Ai = 3.25 and 
A 2 = -4. 

For a slightly different driving frequency /o = 1.879 MHz, we observed a different scenario, 
shown in Fig. 4. In this case a pair frequencies j\ = 1.242 MHz and f% = 0.638 MHz is 
excited above threshold. The pair is also clearly asymmetric in the amplitudes, as shown 
in Fig. 4 (a) corresponding to a driving amplitude of 8 V. Different from the previous case a 
linear growth of the fundamental mode below threshold is observed in Fig. 4(b), saturating 
at a constant value above threshold (7 V in this case) . As discussed in the previous section, 
this situation is typical of a highly detuned second harmonic, where the fully dispersive 
model represents a good approach. The theoretical bifurcation diagrams in Fig.4(b) were 
obtained for the set of detunings A = 1.35, Ax = 0.6 and A 2 = 10, and show a good 
agreement with the measured values. 



C. Self-modulation 



Further increasing the pump above the threshold the field components often develope 
low frequency side bands, as shown in Fig.5. The parameters are those corresponding to 
Fig. 4, and the pump parameter is nearly twice the parametric threshold value. In the 
temporal domain, the presence of the side bands implies a low frequency modulation of the 
field amplitudes. This phenomenon, also demonstrated in solids (15] is usually called self- 
modulation since the source amplitude is kept constant. Self-modulation can result from 
different mechanisms, as dicussed in 



In 



16|. 



14j we considered the stability of the solutions of Eqs.flH]) and (jUJ), identifying the 



appearance the self- modulation with the Hopf bifurcation that the amplitudes undergo under 
certain conditions. The frequency of the slow modulation predicted by the theory was of 
the order of several kilohertz, for typical operation conditions. This is in agreement with 
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the observed side band frequency, Af = 34 kHz, in Fig. 5. A careful analysis of the 
spectrum allows to identify these new frequencies as the result of the mixing (sum and 
difference, /j ± fj) between the fundamental and the subharmonics components, with the 
second harmonics 2 fa, owing to the quadratic nonlinearity of the medium. Some of these 
frequencies are correspondingly labeled in Fig. 5. The experiment also shows that the 
appearance of the low side band frequencies occurs once th pump reaches a critical value. 
This behaviour of the spectrum is chracteristic of a secondary instability of the Hopf type, 
in agreement with the analysis performed in [jjj]. 



IV. CONCLUSIONS 



The models proposed previously to describe the parametric sound generation in acoustic 
resonators can not account for a number of effects observed in the experiment. The first 
reasonable approach to this problem is to consider the role played by the first higher har- 
monics, whose magnitude is strongly dependent of the distance to the cavity resonances. 
The proposed model, which is a simple extension of previous (dispersive) models, consid- 
ers the coupling of the fundamental field to its second harmonic. We demonstrate that 
the bifurcation diagrams predicted by the theory and observed in the experiment can be 
qualitatively different under different parameters, and are very sensible to the experimental 
conditions (mainly the driving frequency and the temperature, both affecting the detuning 
values). Such model also predicts the emergence of secondary instabilities leading to time 
self-modulated solutions, in agreement with the experimental observation at drivings be- 
yond the threshold. We expect that, at even higher drivings, the dynamics of the system be 
chaotic, as the theory predicts. Experimental work in the direction is in progress. 
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a. Figure Captions Figure 1. Schematic representation of the field and cavity spectra, 
in the cases of (a) ideal (loseless) cavity, (b) a real cavity with finite impedances. 

Figure 2. Development of bistability of the fundamental mode as the second harmonic 
detuning is varied. Parameters are 7 = 1,A = —2, and A 2 = 1 (a), A 2 = (b), and 
A 2 = -1 (c). 

Figure 3. Experimental bifurcation diagram (symbols) obtained for a driving frequency 
fo = 1.991 MHz. Solid and dashed lines correspond to the fundamental and subharmonic 
amplitudes, as obtained from Eqs. (5) and (8) for A = —2.35, A x = 3.25 and A 2 = —4. 

Figure 4. Experimental bifurcation diagram (symbols) obtained for a driving frequency 
/o = 1-879 MHz. Solid and dashed lines correspond to the fundamental and subharmonic 
amplitudes, as obtained from Eqs. (5) and (8) for A = 1.35, A x = 0.6 and A 2 = 10. 

Figure 5. Development of low- frequency side bands for the experimental conditions in 
Fig. 4, for a driving amplitude V in = 17V. 
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